THERMODYNAMIC OPTIMIZATION OF 
FIN GEOMETRY 



By 

SHAMBHU NATH SHUKLA 


1 



I 54" “S. K “Y p 

jyt 

6MU 

yHE- 



DEPARTMENT OF MECHANICAL ENGINEERINO 

INDIAN INSTITUTE OF TECHNOLOpY KANPUR 

AUGUST, 1984 



THERMODYNAMIC OPTIMIZATION OF 

FIN GEOMETRY 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

SHAMBHU NATH SHUKLA 


to the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

AUGUST, 1984 





SL 3^ i 


friE^ ' S-HU-THE 


11 


CERTIFICATE 

This is to certify that the work on Thermodynamic 
Optimization of Fin Geometry," by Shambhu Nath Shukla, has 
been carried out under my supervision and has not been 
submitted elsewhere for a degree. 



August , 1984 



( P. 'n. Kaul ) 

Assistant Professor 
Dept, of Mechanical Engineering 
Indian Institute of Technology 
Kanpur 208016 


At- ... 

r«* „ 

iu 




iii 


ACKNOWLEDGEMENT 


I feel very much indebted to Dr. P.N 
wish to acknowledge gratefully the help and 
given by him throughout. 


Kaul and 
guidance 


Shambhu Nath Shukla 



iv 


CONTENTS 

Chapter 

LIST OF FIGURES v 

NOMENCLATURE vi 

ABSTRACT viii 

I. INTRODUCTION 1 

1.1 Introduction 1 

1.1.1 Exergy Losses and Irreversibility 2 

1.1.2 Exergy Loss in Heat Exchangers 2 

1.1.3 Loss to Surroundings 3 

1.2 Literature Review ^ 

1.3 Scope of Present Work 6 

II. MATHEMATICAL FORMULATION 8 

2.1 Entropy Generation Due to Convective Heat 

Transfer from a Single Fin 8 

2.2 Pin Fin 11 

2.3 Rectangular Fin 17 

2.4 Longitudinal Fin of Triangular Profile 23 

III. RESULTS AND DISCUSSIONS 29 

3.1 Results and Discussions 29 

3.2 Conclusions 32 

3.3 Suggestions for Future Work 34 

REFERENCES 35 


FIGURES 



V 


Figure No. 

1 . 

2 . 

3 . 

4 . 

5 . 

6 . 

7 - 

8 . 

9 , 

10 . 

11 . 

12 . 

13 . 

14 . 

15 . 


LIST OF FIGURES 

Details 

Schematic diagram of a general fin in a forced 
convective heat transfer arrangement 

Pin fin 

Rectangular fin 

Longitudinal fin vdth triangular profile 
Ng vs Re^ for pin fin 

Optimum fin diameter for minimum irreversibility 

Optimum pin fin length for minimum irreversibility 

Optimum fin diameter vs parameter, M 

Optimum pin fin length vs parameter, M 

Effect of frictional irreversibility parameter, B 
on the optimum width of the rectangular fin 

Effect of frictional irreversibility parameter, B 
on the optimum length of the rectangular fin 

Entropy generation number, N , vs slenderness 
ratio, Y for a triangular fin 

Optimum width of triangular fin vs. slenderness 
ratio, 

Optimum width vs ^ for triangular fin 
Optimum width vs. B^ for triangular fin 



NOMENCLATURE 


Frictional irreversibility parameter^ equation (2.2.10) 
Frictional irreversibility parameter, equation (2,4,23). 
Drag coefficient. 

Skin friction coefficient. 

Pin fin diameter. 

Drag force. 

Heat transfer coefficient. 

Thermal conductivity of fluid 
Thermal conductivity of fin material. 

Fin length. 

Pin fin conduction parameter, equation (2.2.2). 
Rectangular fin conduction parameter, equation (2.3.2). 
Triangular fin conduction parameter, equation (2,4.6). 
Parameter, equation (2.2.11). 

Parameter, equation (2.4.21). 

Entropy generation number. 

Nusselt number. 

Property group, 

Prandtl number. 

Base heat transfer rate. 

Heat transfer from control volume to remaining fluid. 
Reynolds number, • 

Reynolds number, U^L/ 2 ; . 



Re 


W 


Re 

Re 


s 


gen 

1 

■B 


T 

a 

U 

U 

o 

V 

w 


W. 

X 

y 

a 

K 

C 

S o 

©g 

F 


lost 


Reynolds number, 

Reynolds number, 

Reynolds number, U^S^/ 2>* 

Entropy generation rate. 

Absolute temperature of fin base. 

Absolute temperature of environment. 

1 

Absolute temperature of free stream. 
Parameter, equation (2.4.13). 

Free stream velocity. 

Parameter, equation (2.3.21). 

Width of fin. 

Rate of lost available work. 
Longitudinal coordinate . 

Parameter, equation (2.3.14). 

Half angle, Fig. 4. 

Slenderness ratio. 

Thickness of rectangular fin. 

Base thickness of triangular fin. 
Base-stream temperature difference, T 
Kinematic viscosity of fluid. 

Density of fluid. 



viii 


ABSTRACT 

A theoretical framework for the minimiJm entropy 
generation in extended surfaces or fins has been developed. 
The entropy generation equation comprises two tenns, one due 
to heat transfer and the other due to fluid friction. The 
optimum geometry for a fin is obtained when the above two 
terms add to yield a minimum. Analytical solutions for 
three types of fins viz., a pin fin, a rectangular fin and 
a fin with a trangular profile are presented herein and the 
results plotted on graphs to determine their optimum dimen- 
sions. The thermodynamic optimization has been carried out 
for a given heat transfer rate and flow velocity. 



CHAPTER I 


INTRODUCTION 


1.1 INTRODUCTIONS 

In the design of a system involving the generation or 
use of energy, the exergy method will provide the information 
to better select the component designs and their operating 
procedures that will be most effective for energy conservation. 
The minimization of lost work will produce the most efficient 
system. The exergy method of analysis is based on this princi- 
ple. In this method, the primary task is to determine where 
the losses in the available work occur in the system, what the 
magnitude of these losses is, and which of the losses, when 
corrected, can most effectively improve the system efficiency. 
The selection of the losses in the system to be corrected 
should usually be based on economic considerations. However, 
the use of the avail able -work principle in the exergy method 
will bring about a significant improvement in the data used for 
making judgements, since it takes into account both the quantity 
and the quality of the inefficient energy utilization [l]. 
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1.1.1 Exergy Losses and Irreversibility i 

The exergy losses are produced in a system hy the irre- 
versible production of entropy caused by the non-ideal perfor- 
mance inherent in all real systems and components. Minimization 
of these irreversibilities is necessary in those places where 
the exergy loss to the system is quite significant. This aspect 
of the exergy-analysis method requires a knowledge of the degree 
of irreversibility and the associated exergy losses with mass 
and heat leakage, fluid-flow turbulence and friction. The 
reason for considering all exergy losses in a system during 
initial stages is that some processes that usually have a low 
exergy loss may show upto be the cause of a significant loss in 
another type of system. All irreversibilities in processes 
and systems should be given due consideration if effective 
energy-conservation measures are to be adopted. 

1.1.2 Exergy Loss in Heat Exchangers ^ 

Heat exchangers serve as components in a wide range of 
power and refrigeration equipment. In order to conserve avail- 
able work, it is necessary to approach the design of heat 
exchangers from the point of view of entropy minimization. It 
is important to go further and focus attention on the components 
of heat exchangers and consider the design of each such building 
block for minimum irreversibility. Most of the heat exchangers 
used to-day are generally inefficient from an exergy stand point 
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because they are primarily designed on the basis of low cost 
that dictates a minimum-sized unit. To achieve this small 
sized heat exchanger goal, the temperature difference between 
the fluid streams is maximized. However, the larger the tem- 
perature difference in a heat exchanger, the greater is the 
exergy loss during heat transfer. Hence from exergy- saving 
point of view such designs are not desirable. 


1 . 1.3 


Loss to Surrounding ss 


The loss of thermal energy from a system to the surround- 
ing environment is a direct loss of the exergy in the system. 
This loss is considered in the exergy method of system analysis 
as an internal loss and is, therefore, included in system 
efficiency calculations. 


The first and second laws of thermodynamics, taken to- 
gether, state that the entropy generated by any engineering 
system is proportional to the work lost irreversibly by the 
system [l]. This is expressed concisely as 


W. 


lost 


T Z S 
o gen 


( 1 . 1 . 3 . 1 ) 


where, ^^^ost lost available work or unavailability, 

T^ is the absolute temperature of the environment, and 
Sgen is the entropy generated in each component of the 
system. 
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Equation (l. 1.3.1) implies that the thermodynamic 
irreversibility (entropy generation) of each system component 
contributes to the aggregate loss of available work in the 
system Therefore, it is important to concentrate on 

each component of the system and, by design, try to minimize 
the irreversibility (S ) of that component, 

gGxl 


1.2 LITERATURE REVIEWi 

The traditional approach to the optimization of fins 
consists of minimizing the consumption of fin material for 
the execution of a specified task, Duffin [2], on the basis of 
variational calculus proved that a two-dimensional fin must 
have a parabolic profile if it is to require the least material 
(volume) for a certain heat transfer rate. This design princi- 
ple has been brought closer to reality by Razelos and Imre L3] 
who analysed the role of radiation, two-dimensional heat transfer, 
temperature dependent thermal conductivity and variable heat 
transfer coefficient. Many of these contributions have been 
summarised by Kern and Kraus [4], The essence and practical 
limitations of this design philosophy have been discussed by 
Kraus and Snider [5]. 

An entirely different approach to the design of heat 
exchangers is based on the principle of minimization of exergy 
loss by the system sub-components. This approach was adopted 
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by A. Bejan [6] in the design of counter flow heat-exchangers 
for gas to gas application. His design is based on the thermo- 
dynamic irreversibility rate or useful power no longer available 
as a result of heat exchanger frictional pressure drop and 
stream to stream temperature differences. The irreversibility 
concept is used to establish a direct relationship between the 
heat exchanger design parameters and the useful power wasted 
due to heat exchanger non-ideality. Bejan's paper presents a 
heat exchanger design for a fixed irreversibility. In contrast 
to the traditional design procedures, the amount of the heat 
transferred between the streams and the pumping power for each 
side become the outputs of this design approach, D, Potolikakos 
and A. Bejan [lO] found the fin geometry for minimimi entropy 
generation in forced convection. They derived a general expres- 
sion for the entropy generation rate for an arbitrary fin and 
minimizing this, obtained fin geometric parameters for various 
shapes corresponding to the minimum entropy generation rate. 

In their analysis, they assumed fin-base stream temperature 
difference Qg to be very small con^ared to the absolute value 
of the free stream temperature. This assumption simplifies the 
expression for the entropy generation rate. In the present 
study, however, this assumption has been relaxed and an exact 
expression for the entropy generation rate has been used to 
optimize the fin geometry. 
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1,3 SCOPE OF PRESENT WORK? 

This thesis outlines the procedure for the design of 
fins for minimum entropy generation in forced convection heat 
transfer. This design philosophy allows to properly account 
for the fact that fins, in addition to enhancing the heat 
transfer, increase the fl lid friction also. The irreversi- 
bility minimization philosophy places this trade-off between 
heat transfer and fluid friction on solid foundations as heat 
transfer and fluid friction drag are both mechanisms for entropy 
generation. The competition between these two mechanisms, i.e., 
enhanced thermal contact and fluid friction is settled when the 
heat transfer irreversibility and the fluid friction irrever- 
sibility add to yield a minimum rate of entropy generation for 
the fin. 

The thesis presents a derivation of the formula for the 
rate of the entropy generation for an arbitrary fin engaged in 
forced convection heat transfer. Based on the results therefrom 
it is shown how the geometric parameters of common fin shapes 
can be selected so that the fin saves most of the exergy while 
performing its specified heat transfer function. Throughout 
this study, the classical fin heat transfer model [4j has 
been adopted whereby the fin is slender enough so that the 
conduction process can be regarded as unidirectional. It is 
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further assumed that the properties of the fin material and 
those of the external fluid are constant. The external flow 
is assumed uniform and parallel to the base of the fin. 


CHAPTER II 


MATHEMATICAL FORJIULATION 

2.1 ENTROPY GENERATION DUE TO CONVECTIVE HEAT TRANSFER 
FROM A SINGLE FINs 

The entropy generated hy a single fin in the cross flow 
can be evaluated based on the general model presented in Fig. 1. 
Considering, an arbitrary fin suspended in a uniform stream 
with the velocity, U^, and temperature, T^, the heat transfer- 
rate, qg, is caused by the temperature difference between the 
fin base, Tg, and the free stream,!^. In addition, the cross- 
flow arrangement is responsible for a net drag force, which 
is transmitted through the fin to the base wall. 

As shown in Fig, 1, we choose a control volume which is 
fixed relative to the fluid. The environment (the wall) drags 
the fin with a speed, U^, through the stationary fluid, by 
applying a force % tangentially to the control surface. In the 
steady state, the first and second laws of thermodynamics 
dictate, respectively, 

1B - 'U Fd = 0 

„ "I* % 

Sgen = fZ-’f; 


> 


0 


( 2 . 1 . 1 ) 

( 2 . 1 . 2 ) 
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where is the net heat transfer rate from the control volume 

to the rest of the flioid, and S _ is the rate of entropy 

’ gen 

generation associated with the heat and fluid flow arrangement, 
qg is the heat flow rate at the base of fin and represents 

heat dissipation due to friction. 


Combining equations (2.1.1) and (2.1,2) yields 




gen 


^B 


3 


gen 


cx> 

1 

1 

% ( ^ 
oo 

- i- ) . 

Tb ’ X- 


Fn U 

D 


T 

OO 

i “ Too) 

is the temperature ' 

fin and the free stream. 

(2.1.3) 

can be rewritten as 

% 



(2,1.3) 


(2,1.4) 




The entropy generation rate, given by equation (2,1,4) 
demonstrates that inadequate thermal condi^tance and fluid 
friction contribute hand-in-hand to the degradation of the fin 
thermodynamic performance. The first term on the right hand 
side of the equation represents the entropy generation due 
to heat transfer across a non-zero temperature difference. 
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while the second term is the entropy generated due to fluid 
friction. 


^gen “ ^ gen ^ he at 


transfer 


^^gen"^ fluid 


friction 


(2.1.5, 


In order to minimize the heat transfer contribution to 

S , one must minimize the base stream temperature difference 

practical terms, however, to minimise Qg would imply 

the use of an infinitely large fin for a constant heat flow 

rate, such fin would be impractical and thermodynamically 

undesirable because an infinitely large fin would have an 

infinite (S ) fluid friction, hence an infinite total S . 
gen 

Thus, although, fin size influences monotonically the base 
stream temperature difference, it plays a crucial trade-off 
role in the thermodynamic performance of the fin. 

In this thesis, the fin size trade-off by minimizing 
the rate of entropy generation has been examined in detail for 
some of the most common fin geometries encountered in practice. 

The size trade-offs presented later are based on minimiz- 
ing equation (2.1.4) subject to following constraintss 

1. Constant q^, T^, , 

2. Constant fluid properties. 

Constant fin material properties. 


5 . 
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2.2 PIN FIN 


Pin fin geometry depends on only two dimensions i the 
length, L, and the diameter of the circular cross-section, D. 
According to unidirectionel heat conduction model described in 
the introduction the relationship between the base heat flow 
rate and the base stream t^:mperature difference is [4] 


0 


B 


^ kg ffi tanh (mL) 


( 2 . 2 . 1 ) 


where , 

m * (4h/k (2.2.2) 

s 

and h and afe the heat transfer coefficient and the thermal 
conductivity of the fin material, respectively. 


Substituting this expression for ©g into equation (2.1,4) 
gives the total entropy generation rate as 


\en * \ (niL)(l 


This can be simplified to 


^ T k„D^m tanh (mL) 

4). Od g 


)3 




(2.2.3) 
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Heat transfer coefficient h can be expressed as 

h = (Nu k^/D) = (Nu U^Re^-i;) (2.2.4) 

where is the kinematic viscosity of the fluid, k^ is the 
thermal conductivity of the fluid and Re^ is the Reynolds 
number based on the diameter of the pin fin. 

Now equation (2.2.2) can be simplified to 
m = ( 4 h/kgD )^/2 ^ Uj/(ReQD.kg Rej^O)))^/^ 


Thus 

m = 2(kj/kg)^/^ Nu^/ 2. (U^(Rej 3 ll)) 

(2.2.5) 

and 

L = (Rej^^) /Uj 

(2.2.6) 

Rsj^ is the Reynolds number based on length of the 

fin. 

Force 

Fj^ on the fin is given by- . • 



= Cd-' 

(2.2.7) 

where 

Cq is the drag coefficient. 



Combining equations (2.2.3) to (2.2.7), we get the 
following equation for the entropy generation rate 


\en = T^2(k^kg)^/2N^l/2 ^ tanh(2Nu^/^ 

(k/kg)l/2 Rej^/Rej^) 4- qgT^) Re^ 


( 2 . 2 . 8 ) 



13 


This equation is non-dimensionalised to yield the entropy 
generation number [ 9 ] as 

= (|- (k^/kg)^/^ Nu^/^ Re^ tanh [2Nu^/2(k^/kg)^/^ 

Re^/Re^] + M)“^ + | Re^ Re^ (2.2.9) 

where B is a fixed dimensionless parameter that accounts for 
the importance of the fluid friction irreversibility relative 
to the heat transfer irreversibility, 

B = fi? kg Tyqg2 (2.2.10) 

M is also a dimensionless parameter given by 

- M = (% (2.2.11) 

Parameters B and M are known as soon as fluid properties, 
flow velocity, its temperature and base heat flow rate are 
specified. 

If the fin is slender, the Nusselt humber and the drag 
coefficient can be evaluated from the following results developed 
for a single cylinder in cross flow [ll] 

4 < Re^ < 40, Nu = 0.919 Re^^*^®^ pj,l/3 

Cq = 5.484 ReQ“°*^^^ 
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40 < Rej^ < 4 X 10^, Nu = 0.683 Pr^/3 

Cq = 5.484 Re^“°*^^^ 

4 X 10^ < Rej) < 4 X 10^, Nu = 0.195 Pr^/^ 

Cj^ = 1.1 

4 X 10^ < Re^ < 2 X 10^, Nu = 0.0268 

Cj^ = 1.1 (2.2.12) 

The entropy generation number, N is a function of six 

O 

dimensionless groups, two pertaining to the fin geometry 
(Rej^, Rej^), and four accomting for the working fluid and the 
fin-stream convective arrangement (Pr, kg/k^, B, M) . 

Minimization of N_, with respect to Re. is achieved by 

S ij 

solving 

4Ng/ ^Rej^ = 0. 

Partial differentiation of equation (2.2,9) "with respect 
to Rej^ set equal to zero, yields 

-u(^) Nu Sech^ [2 Nu^/2 (^)^/^ 

^s s D 

+ BCj^ Re^E I (^)^^^ Nu^/^ Re^ tanh [2Nu^/^ 
(^)^/^ + m]2 = 0 

kg Ke^ 

(mL) + B’ tanh (mL) + C' = 




A* tanh^ 


or 


0 


(2.2.13) 
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where, 

A’ = ^ Nu(l + I B Cq Re^^) 

B’ » B Cj^ Re^^ M Nu^/^ 

Ic 

C = I B Cq Re^ - TX (j^) Nu 

s 

,L = 2NuV 2 (^)V2 |i 

Positive solution of equation (2,2.13) is given by 

-B' + (B'^ - 4 A'C 

tanh (mL) = — — (2.2,l4) 

2a' 

Substituting escpressions for A', B’and C’ in (2.2,14) yields 
tanh (mL) = (• B M + [471; ^ Nu(l + — B C^j 

3 

, o T /o Ttk JIu 1/2 

- ^ 3 Re^ ) ( 1 + 

iBC^Rej^^)) (2.2.15) 

From equation (2,2.15), expression for Re opt can be obtained 
as follows 

p— 

mL = tanh”^ [(- '^- B Cj^ Re^^M + (47t ^ Nu (l + 

3 

I B C33 Re^^) - I- B Re^ ^ Nu)^/2(^^ 

s 

|B Cq Rej,5))] 
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or. 


Re 


L,opt 


[ % Nu (1 + I B Cq Re^^) - I B Cj 3 Re^^ 




^ Nu)^/2 (1 ., I b Cq Ren^))] 


'D “'^D 

(2.2.16) 


The engineering significance of equation (2.2.16) is that the 
optimum fin length can be calculated immediately, provided 
Rej^ is specified. Substituting equation (2.2.16) into equation 
(2.2.9)» the entropy generation number N , corresponding to 
optimum pin length is obtained. This new function is dependent 
only on the diameter of the fin as far as the geometric para- 
meters are concerned, and it is minimized numerically as shown 
in Fig. 5, The expression for N for optimum pin length is given 
as follows 

k. 


N 


= [| Nu^/ 2 Rejj [(-'— B Re^^^ M + 




ih% ^ Nu (1+ I B Re^^) - I B Cq Re^ 
s 

(7t ^ Nu)^/^(1 + I B Cj3 Rejj^))] +M] 


'1' 4. I 


B Re 
(Nu) 


W c(- 


Cjj Rejj^M + 


[% Nu(l + I B Cj3 Re/) - I B Cjj Rejj 


(n ^ Nu)^/ 2 I- B Cn Re/))] 


k 


8 


'D 


(2.2.17) 
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Equation (2.2.17) is numerically minimized to give 

^®D,opt substituted back into equation (2.2.16) to 

get Rot , . 

® L,opt 


2.3 RECTANGULAR FINs 


This section focuses attention on fin geometries modelled 
as thin conducting plates parallel to the flow direction. The 
fin geometry analysed here is shown in Fig. 3. The minimimi 
irreversibility design criterion of this fin requires the 
selection of three geometric parameters, viz., the length L, 
the breadth (length swept by fluid), W, and the plate thickness 
S * It is assumed that the length and the breadth of the fin 
are large compared to its thickness, ^ . The relationship 
between the heat transfer rate and the base fluid temperature 
difference is given as [4] 


0. 




B 


k, W m’ tanh (m' L) 

o 


where m' = 


5 


and h is the heat transfer coefficient 

2Nu (k^/kg) 


m’L = [ 


Rog Rg-^ 


,]l/2 


Re, 


where 


Nu = ~ ; Re^ 

^f b 


U 


U 


j Re 


W 


U W 

oo 


(2.3.1) 

( 2 . 3 . 2 ) 


(2.3.3) 


U L 

and Re. = 
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The laminar heat transfer and the skin friction results 
for two-dimensional flat plates as given in [l2] are 

Nu = = o:664 Pr^/^ (2.3.4) 

C„ = 1.328 Re (2.3.5) 

Using equations (2.1.4) and (2.3.1) the entropy gene- 
ration rate for a rectangular fin in the laminar flow region 
becomes 

®gen " ^ 1 b ^oo) 

( 2 . 3 . 6 ) 

The entropy generation number is given by 

N = S /(q^% /k ^>T ^) 
s gen s °° 

= [(2Nu Re^ Re^p^/^ tanh (m'b) + M]"^ 

+ BC^ Re^ Rej^ (2.3.7) 

Combination of equations (2.3.3) and (2,3.4) gives 

m'L = 1.15 (^)^/^ Pr^/^ Re^ (2.3.8) 

s ° 

It is seen that the thickness, g , appears only in the 

heat transfer term of the K equation (2.3.7), consequently, 

O 

Rq^ does not play any role in the minimization process of N^. 
Also, since in most practi'-al cases, £ is determined by 
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consideration such as price, availability and machinability 
of sheet metal, it makes engineering sense to regard Re^ as 
fixed. The minimization of with respect to Re.^- and Re^ 
is achieved by solving the follov/ing simultaneous set of 
equations 


= ° 

Differentiating (2.3.7) 
equal to zero, yields 


and 


tRei. 


with respect to Re^ 


0 

and setting it 


a” tanh^ (m'L) + B" tanh (m’L)+C'' = 0 (2.3.9) 


where, 

A" = 2Nu ^ (1 + B Reg Re/) 

B" = 2 2 B M (Nu Reg )^/^ 

and C " = B Re,,^ _ 2 ~ Nu 

3 

Solution of the above equation (2.3.9) is given by 
tanh (m’L) = 


( B"^ - 4A"C'’)^^^ 
2A " 


( 2 . 3 . 10 ) 


Substituting e 3 ?pressions for A", B” and C" we get the following 
expression 


tanh (m'L) = [-BC^ ^ Nu(l + 

s 

BC^ Re^ Re/) - BC^ Re^2]^/2]/( (2Nu 

(1+BCf R^ Be/)) 
s 


(2.3.11) 
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Combining equation ( 2 . 3 . 8 ) with equation (2.3-5) yields 


Re 


L, opt 


= ((Re^Re^)/(2Nu tanh"^((-BC^ Re^^/^ 


k 

k 


f 


Re^^^/^ M + (2 ^ Nu(l + BC^ Re^ Re,/) - 

3 

BC^ Re^^ m2)1/2)/((2Nu “^)^/^ (l + BC^ Re Re/) 

s ^ 

( 2 . 3 . 12 ) 

Substituting expressions for Nu and from (2.3.4) and (2,3.5) 9 
respectively 5 further simplifies the equation ( 2 . 3 . 12 )^-and we 
get 

-1 / ~i 

(2.3.13) 

J 

“f 

where y = tanh (ra'L) 


Re ( tanh-1 (y) 

B,opt 2 Nu k /k, 


= (-1.15 BI'^ Re//^ [(1+1.328 B Re^ Re//^)- 

(2.3.14) 


B]y[2p2]l/2)/(i 2_,328 Re^Re//^) 

P = (kg/k^)^/VPr^/^ 

Similarly equation (2.3.7) combined with (2.3.4) and (2.3=5) 
yields 


N 


= P/[l.l5 Re^^'^^ tanh (m'L) + P M] 


s "''W 


+ 1.328 B Re//^ Rej^ 


( 2 . 3 . 15 ) 


?>N 


Now, 


W 


0 is solved. 
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Differentiating (2,3.15) with respect to Re^^ and setting it 
equal to zero, we get 


Re 


L,opt 


ll.5(— “5) tanh (m'L)]/[^ 


sech^ (m'L) + B[l.l5 Re^^/^ tanh(m 'L) 

o 


+ P Mj^] 


( 2 , 3 . 16 ) 


Substitution of th''- value of Nu from equation (2,3*^) 
into the equation (2.3.13), on simplification gives 


Re, 


'L, opt 

Equating (2.3.16) and (2.3.17), we get 


Re ^1'* Re 1 (l+i', 

1.15 w 2 ^i-y^ 


rrfs I- 1" (2.3.17) 


= (1-5(57^5) Re^^/2 ReJ-/^ tanh (m^L))/(( ^ sech^ (m'L) 


'W 


B [1.15 ^ (2,3.18) 


From (2,3.14) 


tanh (m'L) = y 

p p 

and sech (m'L) = 1 - y 

Therefore, equation (2.3.18) further simplifies to 

ln((l+y)/(l-y)) = 3'y/(|- (l-y^) + B(l.l5 R^^Re^^/\+ PM 




fh 


(2.3.19) 


.imtM 
■“TTrTi^ 

i' ® # •! J ^ 


hm, m. 
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For given values of B, M, P and » equation (2,3*19) 


yields the optimum value of Re^j i.e. Re-, 


Equation (2*3*19) 


'W '^Wjopt* 

has been solved using an iterative technique. An initial value 
for Re^ is assumed and is substituted into (2,3*1'^) to get y. 
This value of y is put into (2.3*19) to see if this equation 
is satisfied. The value of Re^ that satisfies this equation 
is Re,^ opt* opt known, it is then substituted into 

(2.3.12) to obtain the value of Re^ opt* 

A special case arises when M becomes zero. The equation 
(2,3*17) then simplifies to 

y = tanh (m'L) = (v/(l+v))^^^ 


where, v = (1.328 Re^ Re-^^'^^)”^ 


( 2 . 3 * 20 ) 

( 2 . 3 * 21 ) 


Equation (2.3*19) also simplifies to 


ln((l+y)/(l-y) ) = 3y/(0.5(l-y^) + y’^/v) (2.3*22) 

Combining equations (2.3*20) and (2.3*22) leads to an implicit 
result of the form 


In ((v+l)^/"^ + v^/^) = (v + 1) 


. 1/2 


or 


r ^ _I v^/^ + (v+l)^/^) -|2 

^ ^ “ (v.l)^ 


1/2 


(2.3.23) 


Numerical solution to equation (2.3.23) is v = 0,00905. 

Using this value of v in equation (2,3.21), we get 


^W,opt = 18*5 


-2/3 


(2.3.24) 
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Substituting the value of v in (2.3*20) and combining 
it with (2.3.13) and (2.3*24) yields 

Re, , = 0.175 P Re.^/^ (2.3*25) 

L, opt ' S 

The initial value assumed for Re-^ for the solution of 

equation (2.3.19) was obtained from equation (2.3*24). Figs. 10 

and 11 summarize the numerical results for the optimum width 

of the fin Re,^ and its optimum length Re^^ opt* 

2.4 LONGITUDINAL FIN OF TRIMGULAR PROFILE? 


Triangular fin has practical importance because it gives 
maximum heat flow per unit weight with ease of manufacture. For 
a longitudinal fin with arbitrary profile, the governing diffe- 
rential equation is [4] 


fn (x) ““Y 
dx 


d^O ^1^^) dO 2h 


dx dx 


k. 


Q 


= 0 


(2.4.1) 


where f2 ^(x) is the cross-sectional area of the fin, 0 is the 
temperature difference between a point on the fin surface and 
the surroundings, h and k^ are the surface heat transfer 
coefficient and the thermal conductivity of the fin material, 
respectively. For a unit width of the fin, the cross-sectional 
area function f^(x) may be replaced by a profile curve function 
f2(x) such that for a unit width 

f^(x) = 2 f 2 (x). 
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The general differential equation then becomes 


-F ^ ^ ^ 

f^Cx) ■ — 5- + 


df2(x) 

___ 


d9 

dx 


JL 

k 


© = 0 


(2.4,2) 


For a triangular fin, the profile function can be 
written as 


f2(x) 


bo 


(x/L) 


(2.4.3) 


where is the thickness of the fin at the base and L is the 
length of the fin. Tip of the fin is taken as the origin. 
Differentiating (2.4.3) with respect to x, we get 


dx 2L 


(2.4,4) 


Substituting this into the governing differential equation for 
the temperature excess 6, equation (2.4.2) yields 

x~ + ||-m’' 2 L© = 0 (2.4.5) 

dx*^ 

where m" = (2h/k 6^,)^^^ (2.4.6) 

o ^ 


The general solution of the above equation is 

0 = (2m"(xL)^/^) + C 2 K^ ( 2 m ■’(xL)^/^) 

where and ©2 are the arbitrary constants and T^(2m "(xL)V 2) 
and K^(2m "(xL)^/^) are the modified Bessel functions of the 
first and second kind of the zeroth order, respectively. 
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In order to have a finite temperature excess 0, at the 
tip of the fin where x = 0, the arbitrary constant C 2 must 
equal zero because becomes unbounded. This leaves 

0 = (2m"(xL)^/^) ( 2 . 4 . 7 ) 

is next evaluated using the following boundary condition at 
the base of the fin i.e. 


at X = L, 0 = 0g 


(2.4.8) 


The solution of the equation (2,4.7) is 


I (2m "(xL)^/^) 

0 =: 0p, n. . II.— mr . iii, ■ ii 

Iq (2m"L) 

Heat flow rate for a unit width of the fin is given by 


( 2 . 4 . 9 ) 


^ I 

s dx ' x=L 


(2.4.10) 


Evaluating the derivative of 0 at x = L and substituting it 
into the equation (2.4.10) gives 


= 2h Og {2m"L)/im'\ (2m "L)) (2.4.11) 

For a width W of the fin, the heat flow rate q^ becomes 

Qg = 2h 0g W I^(u)/(m 'T^^(u)) (2.4.12) 

where u = 2m"L (2.4.15) 

Gg = qg m"lQ(u)/(2h ¥ (u)) 


(2.4.14) 
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The minimum entropy generation analysis for a triangular 
fin can be done on lines similar to those adopted for the rec- 
tangular fin. However, the triangular fins involve an addi- 
tional geometric parameter, namely, the half angle, a. 


Assuming a laminar boundary layer flow, the follov/ing 
correlations have been used 


and 


p 'U ^ l;/ 

r = c 

D f cos a 

Nu k 


= 1.328 Re^.^ 


k, 


- 1/2 


h 


= 0.664 Re. PrV3 ^ 


(2.4.15) 


(2.4.16) 


Combining (2.4,6) with (2.4.16) gives 
u = 2m "L = (2h/k^ L 


s '^O' 


where Re 


2.30 (k^/kg)^^^ Pr^/^ Re^^ (2.4.17) 

UcoSc 


w 


u w 


U L 

Re o 


CO Oq 


Substitution of (2.4.14) into equation (2.1.4) yields 

s^... = + Rb T„) + Fj, U^/T„ (2.4.18) 


°o m 


gen ^ -Q 

From equations (2.4.6j and (2.4.16) 


2hW 


1.15(k^ 


In non-dimensional form, the entropy generation number can 
be expressed as 
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= [1.15 Prl/6 Re/'* Re,^/^ * I> 

I'J-I 


(2.4.20) 

where M' = 

(2.4.21) 


Substituting for in the above equation from equation (2.4,15) » 
we get 


N 


(1.15 Pr^/6 ^ ^ 

+ 1.328 Re^^/cos a (2.4.22) 


where 


B. 


1 


= T„ (k^ 


(2.4.23) 


If a is the half angle of the fin and Y ^the slenderness ratio 

i . e. 


y = L/¥ (2.4.24) 

then from the fin geometry (Fig. 4) 

Re^ = yRe^^ (2.4.25) 

and Re ^ = 2y Re^^ tan a (2.4,26) 

Making these substitutions into equations (2.4.17, 2.4.22), 
we get 

u = 2m "L = 2.30 (Y/2tan a)^/^ Re^^/V? (2.4.27) 
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Ng = (1.15 {2^ ten m' )“^ + 

-LqV ; ^ 

1.328 Re^^/^coi a (2.4.28) 

For the evaluation of the entropy generation number, 
the value of P has been tahen as equal to 100 and that of 
Prandtl number as 0,696. For given values of a, M' and 
Y equation (2.4,28) has been numerically minimized to deter- 
mine the optimum width of the fin i.e. Re,^ 


CHAPTER III 


RESULTS AM) DISCUSSIONS 


Representative results of the numerical work of minimiz- 
ing entropy generation rates and variations in optimum fin 
dimensions with respect to various parameters have been pre- 
sented graphically from Fig. 5 through to Fig. 15. 


Fig. 5 shows the variation in entropy generation 
number (N ) with the diameter (Re^) of the pin fin for constant 
B values. Entropy generation number has a clear minimum with 
respect to fin diameter. Reynolds number, based on the diameter 
of the fin, corresponding to minimum is also 

observed from the curves that as B decreases, the shifts 

towards the right thereby increasing its magnitude. Since B is 
the ratio of irreversibility caused by friction to irreversi- 
bility caused by heat transfer, any decrease in its value means 
a higher rate of heat transfer, q^. Higher q^, corresponding 
to higher ©g, naturally augments the entropy generation rate. 
Increase in as a result of decreasing B reduces the 

value of Og, thereby decreasing the total entropy generation 

rate. 


Fig. 6 gives the change in with B for different 

values of M. As discussed earlier, with increasing B, 
decreases. However, the effect of M does not look significant 
of its values considered here. 


for the range 
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Fig. 7 presents the change in the optimum pin fin length, 

calculated by substituting into equation (2.2.16), with 

B for various M values. Re^^^p^ is also seen to decrease with 

increasing B. However, M does not seem to have any significant 

effect on Rot 

L,opt 

The information presented in Figs. 6 and 7 has been 
redraivn in Figs. 8 and 9 with a view to bring out the effect 
of parameter M on and Re^^ opt’ respectively. These 

figures show clearly that variation in M does not influence 
much the optimum values of diameter and length of the fin for 
fixed values of B, However, the slight decrease observed in the 
values of Re^^^p^ and in the above figures towards 

higher values of M can be explained as follows. Since B is 
fixed for a particular combination of fin metal, fluid and 
heat flow rate, q^, it follows that an increase in M can result 
only in enhancing the fluid velocity. Increased flow velocity 
not only increases the heat transfer rate but also the fric- 
tional irreversibility of the system. To maintain a desired 
heat transfer rate and to reduce the frictional irreversibility 
under increased flow velocity conditions, the optimum fin 
dimensions tend to decrease. 
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Fig. 10 shows the variation in optimum fin width for 
a rectangular fin with frictional irreversibility parameter, 

B, for constant M, It is found that the optimimi fin width 
decreases as frictional irreversibility parameter, B, becomes 
more pronounced. -Again for a given value of B, optimum width 
of the fin decreases vd.th increasing M, 

Fig. 11 presents the variation of optimum fin length 
of rectangular fin with frictional irreversibility parameter, 
B, Optimum fin length is seen to decrease with B, Similarly 
for a fixed value of B, there is a decrease in the optimum fin 
length with increase in M but in this case the effect of M is 
much more pronounced. It is to be noted that trends observed 
in the case of rectangular fins are in general agreement with 
the conclusions reached in the case of the pin fin. 

Fig. 12 gives the variation of with Re^ for a trian- 
gular-profile fin, for various slenderness ratios l/W). 

The fin irreversibility nimiber, N^, assumes a minimum at a 
specific value of fin base width which is designated 

Fig. 13 depicts the variation of various 

slenderness ratios for fixed B^^, values. Regardless of the 
value of B^, it is seen that the optimum width of triangular 
fin decreases as fin becomes more slender. Increase in y 
decreases the base-stream temperature difference, 9g, thereby 
decreasing the entropy generation due to heat transfer. 
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However, from ^pquation (2.4.28) it is clear that entropy gene- 
ration rate due to heat transfer is inversely proportional to 
the square root of ^ whereas the frictional entropy genera- 
tion rate is directly proportional to . The slenderness 
ratio, y , therefore has a more pronounced effect on the fric- 
tional irreversibility. Re-^ therefore, tends to decrease 

with slenderness ratio, y , so that frictional irreversibility 
is reduced. 

Fig. 14 exhibits the effect of M ' on the optimum width 
of the fin. ' For given value of y the optimum width is seen 
to decrease with increasing M’ . This is also in agreement 
with the results obtained in the earlier two cases viz. the 
pin fin and the rectangular fin. 

Fig. 15 delineates the variation of with B for 

various values of a for M’ = 1.0 and ^ = 4,0. It is apparent 
from the curve that the half-angle, a has a relatively small 
impact on the optimum fin size for minimum irreversibility. 
However, the optimiom width of the fin decreases with increase 
in the half -angle, a, because of increased frictional irrever- 
sibility. 

3.2 CONCLUSIONS s 

Main conclusions drawn from this study are as follows 

1, The optimum diameter and the length of the pin fin 
decrease with frictional irreversibility parameter B to 
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reduce the entropy generation rate due to friction, 

2. Parameter M does not have any noticeable effect on 
the optimum geometry of the pin fin, for most of the range 
considered. There is, however, slight decrease in the optimum 
dimensions of the fin towards higher values of M. 

3. The optimum dimensions of the rectangular fin decrease 
with frictional irreversibility parameter, B, In this case, 

M has a more pronounced effect on the optimum length of the 
fin. Both dimensions, i.e., length and width, register a 
decrease with M. 

4. As the triangular fin becomes more slender, its optimum 
width decreases to reduce the effect of friction on the entropy 
generation number. 

5. The optimum width and length of the triangular fin 
decrease with M’, However, this decrease is small, 

6. The half angle a has a relatively minor impact on the 
optimxjm fin size for minimum irreversibility. 
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3.3 SUGGESTIONS FOR FUTURE WORKi 

A similar study can be carried out for a heat exchanger. 
The heat exchanger should be optimized on component-by— component 
basis. 

In the present study, constant fluid and fin material 
properties have been assumed. Thermal dependence of fin 
material and fluid properties e.g., ^ should also 

be considered in the optimization analysis, heat loss from the 
fin tip should be taken into account. 

Thermodynamic optimization could also be studied for 
different types of arrays of fins e.g., regular or staggered 
arrangements. For this experiments will have to be performed 
to determine the looal heat transfer and drag characteristics 
3 ,rray . Also spacing between the fins and other dimtn 
sions will be determined with a view to improve the thermody- 
namic performance of the arrays. 

Similar approach could also be adopted for other shapes 
of fins like trapezoidal fins, fins with parabolic profile, 
conical fins etc. The fin geometric parameters should be 
designed such that minimum irreversioility is produced. 
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